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INTRODUCTION 

There is a relatively stable component of the large-scale corona that persists for periods of the 
order of a solar rotation. This component is captured in the synoptic charts constructed from the 
daily coronagraph observation of the density structures visible in Thomson-scattered light from 
the photosphere (e.g., Fisher and Sime 1984). The Thomson-scattered light is optically thin and 
the coronagraph images correspond to electron densities integrated along the line of sight. A long 
standing need therefore exists for the development of theoretical models of the corona capable of 
relating the observed integrated density to the actual three-dimensional coronal structures. It is 
the magnetic field and electric current in the highly conducting medium of the corona that order 
the density into various structures. Unfortunately, the coronal magnetic field cannot be measured 
to any useful spatial resolution. Only the longitudinal component of the magnetic field in the 
photosphere can be measured routinely with a magnetograph. An obvious approach for a theoret- 
ical model to adopt is to use the measured photospheric magnetic field as an input and extrapo- 
late for both the magnetic field and density in the corona. Until recently, this approach has been 
feasible only with the potential field model (e.g., Newkirk and Altschuler 1970). The disadvan- 
tage, as is well known, is that a potential magnetic field does not interact with the plasma and 
the model predicts no density structure. The idea to use linear force-free magnetic fields with a 
constant a, proposed by Nakagawa (1973), has difficulties of nonuniqueness and infinite total 
magnetic energy when applied to coronal models in spherical geometry (Seehafer 1978). Non- 
constant a force-free magnetic fields and general non-potential magnetic fields in static equili- 
brium have so far proved to be mathematically intractable in three-dimensional geometry. In this 
paper, we describe a special class of magnetostatic equilibria, which are mathematically simple 
and yet sufficiently versatile so as to fit any arbitrary normal magnetic flux prescribed at the pho- 
tosphere. With these solutions, the corona can be modeled with precisely the same mathematical 
procedure as has previously been done with potential fields. The magnetostatic model predicts, in 
addition to the coronal magnetic field, the three-dimensional coronal density which can be com- 
pared with coronagraph observations. This is an ongoing project being conducted at the High 
Altitude Observatory (HAO), using data from the Kitt Peak Magnetograph and several corona- 
graphs operated by HAO. 

THREE-DIMENSIONAL MAGNETOSTATIC MODELS 

The corona is in a state of expansion (Parker 1963, Hundhausen 1972). The expansion is believed 
to be negligible in the low corona, and, as a first approximation, the low corona is taken to be in 
static equilibrium and described by 

^(vxfl)xS-V P -f^-r - 0 , IX] 
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where B , p and p are, respectively, the magnetic field, pressure and density, G is Newton’s grav- 
Rational constant, Mq the solar mass and r the heliocentric distance. To close the equations, 
we may relate the temperature to p and p through the ideal gas law and specify an equation for 
steady energy balance. The balance of energy in the corona plays an important role in determin- 
ing coronal structures. The detailed processes are complicated and our current knowledge about 
them is very incomplete. In any case, to account for these processes would render the magntos- 
tatic problem quite intractable. It therefore seems reasonable, as a first step, to ignore the energy 
equation and seek solutions to equations (1) and (2) with no further constraints. Even within such 
an artificial consideration, the magnetostatic problem is not readily tractable (see the review by 
Low (1985a)). Equations (1) and (2) have not been treated in fully three-dimensional geometry 
except for the cases of the potential magnetic field and the constant-a force-free magnetic field. 

In a recent development, Low (1985b) showed that equations (1) and (2) can be reduced to a 
simpler problem without assuming any special geometry if the electric current density 7 is 
assumed to flow everywhere perpendicular to the solar gravity. In this case, employing the usual 
spherical coordinates, we can express the magnetic field and electric current density as 
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where $ is an arbitrary function of 
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and r . The magnetostatic problem reduces to solving 
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for <£ as a boundary value problem given 4((-— ,r ). The freedom to choose the form of arises 

dr 

from our not imposing an energy equation. From the solution 4>, the magnetic field is given by 
equation (3) and the equilibrium pressure and density are given by 
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where p 0 and p Q are the pressure and density of an arbitrary spherically symmetric hydrostatic 

d<$> 

atmosphere. The familiar case of a potential magnetic field is recovered if we set = -r— , and 

or 

then only the hydrostatic components in equations (6) and (7) remain. In general, equation (5) is 
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nonlinear and suitable numerical methods need to be developed to solve it. If we set 


*-*>^7, [8] 

where 77 is given, equation (7) is linear in $ and standard numerical methods for elliptic equations 
can be employed to solve the boundary value problem. Bogdan and Low (1986) considered 
several cases of given rj for which the boundary value problem is soluble analytically in terms of 
classical functions. In particular if 


„(r)=(l+f) 2 , [9] 

where a is a constant, equation (7) admits the series solution 

Hr M) = E»°°= 1 Em = -» (>U (r + a )— 1 + B mn (r + a )» ) Y n m (0,<f>) , [10] 

where A mn and B mn are constant coefficients and Y™ are the usual spherical harmonics. As to be 
expected, the classical spherical-harmonic expansion for the potential field is obtained if a = 0 . 
In fact, equation (9) shows that ^ -+ 4> in the limit of r — ► oo so that in that limit, the magnetic 
field is potential. In other words, the electric current is confined to the low corona within a 
heliocentric distance of the order of a . For a given non-zero a , equation ( 10 ) can be fitted with a 
prescribed radial magnetic field at an inner spherical surface identified with the solar surface, 
using B r = An outer boundary condition is needed to determine the coefficients A mn and 
B mn uniquely. This outer condition may take the form of requiring the magnetic field to vanish 
as r — ► oo or that the magnetic field becomes purely radial at some outer spherical surface (e.g., 
Newkirk and Altschuler 1970). The procedure follows essentially the same steps as that of the 
classical potential field model. However, the magnetic field is not potential for a non-zero with 
the equilibrium pressure and density is given by equations ( 6 ) and ( 7 ). 

DISCUSSION 

We plan to generate three-dimensional magnetospheric models of the corona based on equations 
(9) and ( 10 ), with an input of the photospheric normal magnetic flux from the Kitt Peak Magne- 
tograph. This is a first opportunity for constructing a three-dimensional magnetostatic corona 
which is based on photospheric magnetograms and can be tested against coronagraph observa- 
tions. By varying the free parameter a , it will be interesting to see if we can achieve a reason- 
able agreement between the density predicted by the model and the observed line-of-sight 
integrated density. The model is based on several artificial assumptions and we do not have 
apriori reasons to expect excellent agreement. The assumption of the electric current flowing 
perpendicular to gravity is not as severe as it may appear. Such a distribution of the electric 
current may result from stretching a magnetic field radially outward, starting from a potential 
state, as demonstrated by the solutions described by equation ( 10 ) for a >0 (Bogdan and Low 
1986). The outward stretch of the magnetic field may result from the buoyancy of magnetic 
structures. The assumption of equation (9) is more restrictive. It requires that the presence of 
electric current can be characterized by a single parameter, namely, a . The illustrative examples 
in Bogdan and Low (1986) look promising, but it remains to be seen from actual modelling with 
data if this is a workable hypothesis. 
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